Theory of interlayer tunneling in bi-layer quantum Hall ferromagnets 
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Spielman et al. have recently observed a large zero-bias peak in the tunnel conductance of a bi-layer 
system in a quantum Hall ferromagnet state. We argue that disorder-induced topological defects in 
the pseudospin order parameter limit the peak size and destroy the predicted Josephson effect. We 
predict that the peak would be split and shifted by an in-plane magnetic field in a way that maps 
the dispersion relation of the ferromagnet 's Goldstone mode. We also predict resonant structures in 
the DC I-V characteristic under bias by an ac electric field. 



Exotic effects induced by inter-layer Coulomb interac- 
tions have made strongly coupled bi-layer quantum Hall 
systems at total Landau level filling factor v = 1 the sub- 
ject of numerous theoretical and experimental studies. 

When the layers are widely separated they behave 
as two weakly coupled v — 1/2 composite fermion met- 
als. However, when the inter-layer distance d is smaller 
than about twice the magnetic length £, the system spon- 
taneously develops interlayer phase coherence and forms 
H a v — 1 quantum Hall state. This broken symme- 
try state may be described as a Bose condensate [[|]|] 
in a bosonic Chcrn-Simons field theory, as an easy-plane 
ferromagnet P,Q,p|J3] in a theory based on a pseudospin 
representation for the layer degree of freedom, or as a 
superfluid excitonic condensate [lCjJll]] in a theory based 
on a single-layer particle-hole transformation. We use the 
pseudospin language below. 

In a recent experiment, Spielman et al. |l2[ observed a 
qualitative change in the voltage dependence of the inter- 
layer tunneling current I(V) upon entering the ordered 
state. For large d/l, bilayer v = 1 systems exhibit a 
pseudogap behavior: the tunneling current is extremely 
small at low bias voltages. This suppression of tunneling 
is attributed to the slow relaxation of charge character- 
istic of the v — 1/2 state in each layer. In the ordered 
state, Spielman et al. discovered a strong and sharp zero 
bias peak in the differential conductance dl/dV. It ap- 
pears plausible that this peak is related to the Josephson 
effect predicted by Wen and Zee § and by Ezawa and 
Iwazaki. || In contrast to the conventional Josephson ef- 
fect, however, no zero-bias supercurrent (infinite tunnel- 
ing conductance) was observed. The peak conductance, 
though enormously enhanced, did not exceed 10 _2 e 2 /h. 

In this Letter we analyze tunneling in the bi-layer 
quantum Hall v = 1 state. We explain how long 
range density inhomogeneities introduce topological de- 
fects (merons) into the SU (2) pseudospin order parame- 
ter. These defects carry both charge and vorticity |p[Jl3| 
and constitute a dissipative environement which turns 



the Josephson effect into a finite tunneling peak whose 
height and width is a measure of the dynamics of the 
topological defects. We predict dependences of the tun- 
neling current on in-plane magnetic field strength B| | , 
bias voltage frequency, and on the homogeneity of the 
2D layers. In particular, we show that a measurement of 
I{V, B\\) would test the main premise of our theory, the 
existence of one low energy excitation mode, and would 
map the dispersion relation of that mode. Finally, we 
analyze the current distribution for a perfectly homoge- 
neous sample. 

The order parameter field of the quantum Hall fer- 
romagnet is a pseudo-spin unit vector rfi. When fluc- 
tuations out of the easy plane are small, it can be 
parametrized by an angle ip and the conjugate 'charge' 
m z : fh — (cos ip, siny, m z ). In the absence of tunnel- 
ing, disorder and topological defects, the long wavelength 
Hamiltonian density of the v = 1 bi-layer state is [0J|,01 
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where no = is the average density. In Hartree-Fock 
theory p s ~ 0.4K and the capacitance T is increased from 
its electrostatic value. Since the momentum density 
conjugate to p isp^ = TiriQm z /2, the Hamiltonian (Q) has 
a single linearly dispersing collective mode with velocity 
u = yj p s /r. This Goldstone mode signals superfluidity 
for in-plane currents which are antisymmetric in the layer 
index. Hm Taking proper account of the significant ex- 
change enhancement of T yields u ~ 0.1e 2 /he. Eq. (|l|) 
is valid only in the limit q — > 0. Away from that limit 
this collective mode has a more complicated dispersion, 
denoted by u> q , which shows roton effects. [^],[| It is this 
dispersion curve that may be extracted from a measure- 
ment of I(V,B\\). 

The inter-layer tunneling operators are 



T± = - d 2 r A(r) 



e ±i<p(r) e ±iQ 



(2) 
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where the ± sign refers to the direction of tunneling, 
Qb = —f^r is a characteristic wave vector introduced 
by the magnetic field Sm (we choose the gauge A\\ = 
xB\\z). The quantity A = g^p-AsAS is proportional to 
the tunneling amplitude and may vary with position due 
to disorder in the tunnel barrier. Here we do not discuss 
this source of disorder since, on its own, it can not de- 
stroy the Josephson effect. As in a Josephson junction, 
the tunneling term in the Hamiltonian is T + + T_ , while 
the tunneling current operator is ie(T + — T—)/h. The 
striking similarity of the expressions above to their coun- 
terparts in superconducting Josephson junctions make it 
clear that a calculation of the tunneling conductance un- 
der Eqs. ([[]) and (||) leads to a Josephson effect, in con- 
trast to the experiment. We now explain the way that 
disorder destroys the Josephson effect in the present sys- 
tem. 

In the v — 1 bi-layer system, a deviation of the to- 
tal density from v = 1 introduces topological defects 
(merons) into the order parameter vector m. In terms 
of bosonic Chern-Simons theory, this statement is a con- 
sequence of the residual magnetic field left, away from 
v = 1, after the external and Hartree-Chern-Simons mag- 
netic fields almost cancel one another. This residual 
field introduces vortices into the bosonic order param- 
eters of the two layers. In the language of a quantum 
Hall ferromagnet, this observation is a consequence [QJ^] 
of the coupling of the symmetric density to the order 
parameter fh. The symmetric part of the density is con- 
strained to satisfy n(r)-no = ■^eabe^K»V<9 a TO l ,<9&m K 
V • (^z x V(f) - x \7ip. The deviation from n is 
then composed of a charge density carried by an electric 
dipole field ^z x Vy, and by a charge density attached 
to topological defects in m. The latter are merons of 
four types, carrying a charge of ±|, and characterized 
by their vorticity (the sign of V x V(f at the core) and 
the layer in which their charge resides (m z at the core). 
Merons interact coulombically due to their charge, and 
by a logarithmic interaction due to their vorticity. Below 
the Kosterlitz-Thouless temperature Tkt ~ Ps, merons 
are bound in pairs of opposite vorticity to avoid the log- 
arithmically diverging energy penalty. 

In realistic samples there are long range density fluc- 
tuations whose relative magnitude is estimated |l4| to 
be 4%. Thus, the typical distance between the disorder- 
induced meron pairs is ~ 12£. The separation between 
the two merons that constitute a pair is estimated to 
be ~ 61, |lj comparable to the spacing among differ- 
ent pairs. This estimate is obtained by balancing the 
Coulomb repulsion and logarithmic attraction. Thus, the 
v = 1 bi-layer sample studied in Ref. is analogous to 
a superconducting junction with random magnetic flux 
that introduces many vortices in the two superconduc- 
tors. Meron pairs may carry a charge ±e (distributed be- 
tween the two layers) or be charge neutral. The charged 



pairs affect the longitudinal resistivity to the flow of sym- 
metric current. In the sample of Spielman et al. this 
resistivity is large (<~ lkfi), indicating that the charged 
vortex pairs are highly mobile. Furthermore, the dissi- 
pation is not frozen out at the lowest attainable temper- 
atures indicating that these objects are disorder- rather 
than thermally-induced. Tunneling in this system is then 
strongly influenced by these merons, in a way discussed 
below. The meron pairs do not however destroy the anti- 
symmetric superfluid mode unless they become unbound. 

Appealing to the experimental observation that there 
is no DC Josephson effect (i.e., current linear in the tun- 
neling amplitude) we may use Fermi's Golden Rule to cal- 
culate the tunneling current perturbatively. For a sample 
of size L 2 

I(V) = 2m ^ L2 [S(Qb, eV) - S(-Q B , -eV)}, (3) 

where S(q,hoj), the spectral density for the fluctu- 
ations of the operator e lv at wavevector q and fre- 
quency u>, is proportional to the Fourier transform of 
^ e jy(r,t) e -i<p(o,o) \ (whgj-g angular brackets denote ther- 
mal average). Our prediction regarding the dependence 
of the tunneling current on B\ \ can now be easily under- 
stood. For weak disorder, the spectral density S(Qb, eV) 
is sharply peaked at 

eV = Thoq b (4) 

Thus, as the parallel field is varied, the peak in the tun- 
neling conductance is shifted in a way that reflects the 
dispersion of the low energy excitation mode. This is 
precisely analogous to the Carlson-Goldman experiment 
measuring the collective oscillations of the pair field in 
a superconductor. |l5f| An observation of this dispersing 
peak would also confirm an essential ingredient of the 
picture we use, namely the existence of a single branch 
of low energy excitations. The parallel field allows only 
tunneling between states that differ by a momentum Qb- 
Energy conservation requires the energy of these states 
to differ by eV. When there is just one low energy excita- 
tion branch, there is only one value of the voltage where 
both these conditions are fulfilled. This is not the case 
for a fermi liquid (for Qb ^ 0). 

To begin our analysis of the effect of merons on 
the tunneling conductance we make the simplifying 
ansatz that the order parameter phase can be sepa- 
rated into the sum a vortex part ip m and an indepen- 
dent spinwave part (p with the former obeying G m (r, t) = 
( e iym(r,t) e -iv> m (o,0)) _ eX p(-^ - The gaussian 

form for the spatial dependence is chosen for algebraic 
convenience. Applying this ansatz to Eqs.(|l|||) yields 

I(V,B n ) = J°°dt J o?rG m {r,t)e-? D ^) 

sm — - — cos Qqx sm — — (5) 
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with a 'Debye- Waller factor' exp(— D(r, t)/2), where 

^(r,t) = ([^(f,t)-^(0,0)] 2 } 

Ehu . . ftug 

[1 - cos(g • r) cos(u^)] coth — (6) 

(we set &b = 1 throughout) and a commutator term 



C(r,t) = i[<p(r,t), ip(0,0)} 
ft 



-0(ttt - r) 



-1/2 



(7) 



which is independent of temperature, and limits the 
r— integral in (|^) to a "light-cone" of r < ut. All cor- 
relators are evaluated in the absence of tunneling. We 
rely on the global £7(1) symmetry and the freedom to 
renormalize £ and t v to partially justify the simplifying 
approximation of neglecting all disorder in the spinwavc 
hamiltonian. 

As long as 27rp s 3> h/r v » T we can expand Eq.(||) 
to first order |l6| in C and approximate D by its zero 
temperature value. In this limit the current becomes 



l(V,B l{ ) = - 



h AT J Pe 



(eV - hw p ) 2 + (<5 P ) 2 (eV + hw p y + {5 v f 



(8) 



where 5 V ee Ti/t v . For large r^,£ Eq. (||) shows a peak 
in the current at the voltage corresponding to the Gold- 
stone mode energy in accordance with Eq.(j|). The effect 
of Tip, £ is to smear this peak over a range of ft/£ in mo- 
mentum and H/t v in voltage. As long as Qb£ 3> 1 and 
uQbt v 3> 1, this smearing is insignificant. 

The expression for the differential conductance sim- 
plifies considerably in the limit Qb — 0, £ <C itr^ and 
eV < * 



dJ _ 1 e z £_ 2 n L 2 A 2 SAS 
dV ~ 8~hP ~p s 6 



{eVy + (8 V ) 2 • 



(9) 



Interestingly, we see that when t v — oo, i.e., when the 
merons provide a random static background phase field, a 
Josephson-like singularity of Mr is still present (as is the 
antisymmetric superfluid property) . As shown below, the 
singularity is present also at finite temperature T <C ps- 
Static topological defects break translational invariance 
and thus open more phase space for excitation of spin 
waves in the tunneling process. However, they do not 
expand the degrees of freedom involved beyond the single 
spin wave mode, and thus do not dephase the process 
enough to destroy the zero-voltage singularity. 

The temperature dependence of (^|) originates from the 
temperature dependence of D and the temperature de- 
pendence of p s and t v . Here we calculate the temperature 
dependence of D. At zero temperature it gives the space 



and time independent result Do = J q £ <s /2 d, 2 ~ ^-8. 
At finite temperature we approximate cothx wl+ \ e ~~ x > 
define dimensionless length and time variables, f = &■ 
and t = Tp, and obtain (for large r, t, and r < ut), 



D{f,t) 
T 



D 



log 



;t + i/2) + V(t + i/2) 5 



(10) 



The temperature dependence of D affects I{V) then only 
at high temperature (T 3> eV), where we can approxi- 
mate t + i/2 k t. For ut v > £ and B\\ = 0, Eq.(|) 
reduces to 



e\ 2 L 2 

Tipsfl 



dt I d rexp( 

n Jr<ut 



n 



'tT\ 



>rT\ 

\ hu ) 



1 

2 U 
eVt 



t 



(11) 



Most of the contribution is then from long times, while 
r is limited to be smaller than £. For a static meron 
background (t v — oo) we find for eV <C T 



di \ 2 eL 2 



pi 



(12) 



which is consistent with the more complete scaling form 
which can be derived in the classical limit from the ex- 
pression of Nelson and Fisher for the dynamical struc- 
ture factor of the XY model. (l7) In the presence of 
a finite t v , the temperature dependence of D affects 
the tunneling current significantly only in the window 
— ,eV < T < Tkt- In the experiment of Spielman et 
si. the peak width is much larger than the temperature. 
Thus, the observed temperature dependence probably 
results from the temperature dependence of p s and t v 
rather than D. 

Using Eq. (||), we can fit the width of the conductance 
peak in the experiment with a phcnomcnological value 
6 V ~ 0.75K. This value gives ut v w 111, which is remark- 
ably close to our estimate of £ based on the meron pair 
spacing. (Sufficiently close that the Lorentzian approxi- 
mation in Eq. (^) for the peak width will be somewhat 
inaccurate.) 

Our naive estimate for the peak height in the experi- 
ment is too large by some two orders of magnitude, but is 
highly uncertain due to the exponential sensitivity to the 
ultraviolet cutoff and the acoustic approximation used in 
computing the Debye- Waller factor. In addition, the es- 
timate As^s « 90/iK is exponentially sensitive to the 
parameters in the modeling of the barrier potential (in 
particular the poorly understood effective mass appropri- 
ate for the high Al concentration in the barrier) and so 
might be off by a significant factor. [Oj It might also be 
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possible that the superfluidity and the tunneling occur 
predominately in isolated regions close to filling factor 
v = 1 containing few vortices. The parasitic series trans- 
port resistance ~ l/a xx in this Corbino-like geometry 
could significantly reduce the peak height. 

We now consider inter-layer tunneling under the com- 
bined effect of a time independent dc voltage and a time 
dependent ac electric field Esinwt, directed perpendic- 
ular to the two layers. As long as the system is not 
heated, this field can be incorporated into our calculation 
by writing T± = - / df \ e ±l ^+^ COSUJt l Repeating 
the calculation carried out above, we find that the tun- 
neling differential conductance 4&(V) exhibits peaks at 
eV = nhuj, with n an integer. This feature is common 
to all tunneling systems where the dc differential conduc- 
tance is strongly peaked around zero voltage (for exam- 
ple, a bi-layer system at zero magnetic field). We note, 
however, that the quantum Hall ferromagnet is relatively 
less prone to heating, due to the small longitudinal con- 
ductivity. 

Finally, we discuss the current distribution in an ideal- 
ized zero-disorder and vortex-free system. Eqs. ([!]) and 
(0) then do indeed lead to a Sine-Gordon equation for 
the phase, as in a long Josephson junction. However, 
due to the two-dimensionality of the problem, the criti- 
cal current is not proportional to the area of the sample. 
Consider a setup where the current is fed into one layer 
from, say, x = — oo, and taken out from the other layer 
at x — oo, and where tunneling is limited to the region 
— 4 < x < -j. Since the symmetric part of the cur- 
rent (I S ym) is conserved, the boundary conditions for the 
Sine-Gordon equation require §f = — ff L=-i = 

I sym . For L > £j = y/4TvPp s /A SA s ~ 4/im the time- 
independent solution to the Sine-Gordon equation in the 

— — \x\ 

tunneling region is f{x) ~ 2 arccos tanh 2 , , tunnel- 
ing takes place only within a distance of order £j of the 
x = ±4 lines, and the maximal current that can tunnel 
is L— independent, and is given by (2e/h)p s W/£ (here W 
is the width of the current contact). For the parameters 
we use this current is ~ 4nA//xm-W ~ 80nA. The ex- 
perimental measurement current was much smaller than 
this value. Thus, the absence of a Josephson effect can- 
not be attributed to a large measurement current. For 
the sample geometry described above and used in the 
experiment, the tunnel resistance is effectively in series 
with the Hall resistance. The observed tunnel resistance 
is however much larger, ~ 10 2 h/e 2 , again indicating that 
there is no Josephson effect. 

To conclude, we have attributed the lack of a Joseph- 
son effect in tunneling measurements in a bi-layer quan- 
tum Hall v = 1 state to density inhomogeneities that 
introduce dynamical topological defects into the order 
parameter. The observed peak width is quantitatively 
consistent with this picture. We showed that a measure- 
ment of the tunneling I{V) dependence on Bn would map 



the dispersion relation of the low energy mode of the sys- 
tem, and that tunneling in the presence of an ac electric 
field would result in resonances at voltages correspond- 
ing to the ac frequency. Finally, we showed that even for 
a perfect sample where the Josephson effect takes place, 
the critical current would not scale with the size of the 
sample. 
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